Generalized Lemaitre-Tolman-Bondi solutions with Pressure 
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Utilizing the ADM equations, we derive a metric and reduced field equations describing a general, 
spherically symmetric perfect fluid. The metric describes both the interior perfect fluid region and 
exterior vacuum Schwarzschild spacetime in a single coordinate patch. The exterior spacetime is 
in generalized Painleve-GuUstrand coordinates which is an infinite class of coordinate systems. In 
the static limit the system reduces to a Tolman-Oppenheimer-Volkoff equation on the interior with 
the exterior in Schwarzschild coordinates. We show the coordinate transformation for the non- 
static cases to comoving coordinates, where the metric is seen to be a direct generalization of the 
Lemaitre-Tolman-Bondi spacetime to include pressures. 
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I. INTRODUCTION 

The Lemaitre-Tolman-Bondi (LTB) class of solutions 
of Einstein's field equations describing spherically sym- 
metric, inhomogeneous dust P, a 13 are well known. 
However, despite many years of work, the generaliza- 
tion to a general perfect fluid^ source is not known. 
There are a number of exact solutions with specific mat- 
ter distributions or geometries, for example the Tolman- 
Oppenheimer-Volkoff (TOV) static models, the FRW 
models with homogeneous pressure and energy density, 
and the list goes on. 

Much work has been done on the collapse of self-similar 
perfect fluids (see for e.g. 0,0 and ;9] for a more recent 
review). The field equations for simple equations of state 
reduce to systems of ordinary differential equations and 
can therefore be analysed in great detail using various 
methods of dynamical systems. Giambo et. al. [Tol [Tlj | 
recently reduced the Einstein field equations for a perfect 
fluid with barotropic equation of state to a single, second 
order quasi-linear differential equation on the metric co- 
efficients. This work was achieved using area-radial co- 
ordinates, and enabled the analysis of naked singularities 
in the collapse of fluids of this type. 

One specific purpose of exact solutions is to describe ei- 
ther a stellar or galactic system. In particular, the equa- 
tions of hydrostatic equilibrium describes the pressure 
required to support these bodies from collapsing under 
their own gravitational field. Furthermore, the non-static 
cases are interesting in that they describe the collapse of 
these systems, and hence the possible formation of black 
holes. 



'Electronic address: paul. lasky@sci. monash. ed u. au| 
tElectronic address: anthony.lun@sci.monash.edu.au' 
^ The definition of perfect fluid we utilize in this article follows 
that of HHH 



However, many known matter solutions only describe 
the region that contains the matter. The exterior re- 
gion is considered as an additional extra, the matching of 
which provides boundary conditions for the collapse pro- 
cess (see for e.g. 0)- However, recently expressed a 
metric that describes both interior and exterior regions 
of a collapsing body as a single solution of the field equa- 
tions. This alleviates the requirement for the matching 
schemes at the interface, a process which is sometimes 
difficult due to the matching of two separate coordinate 
systems. 

In the interior region of the spacetime is a 

marginally bound ball of homogeneous dust, and the 
exterior spacetime is the Schwarzschild spacetime in 
Painleve-Gullstrand coordinates 0, Q|. In we 
generalized this model to non-marginally bound inho- 
mogeneous dust for the interior, with the Scwarzschild 
spacetime in generalized Painleve-Gullstrand coordinate 
system for the exterior region. The coordinate system 
for the dust region of the spacetime is such that shell- 
crossing singularities appear as simple fluid shock waves 

In this paper, we generalize |l5j | to include pressure 
in the form of a perfect fluid. The interior region is first 
derived as a generalization of the coordinates used in |0| , 
and we further show the coordinate transformation to put 
this into comoving coordinates. In this way, the solution 
is shown to be a generalization of the LTB spacetime to 
include pressures. 

By taking the static limit of our coordinates we find 
the system reduces to the TOV equation for hydrostatic 
equilibrium. This solution matches smoothly onto ex- 
terior Schwarzschild coordinates. Considering the min- 
imum radius for a static star is larger than the appar- 
ent horizon, this implies the static system is everywhere 
regular in these coordinates. However, if we consider the 
Schwarzschild coordinates for the collapsing case, then at 
some point through the evolution, the surface of the star 
will be inside the event horizon. Therefore, the spacetime 



2 



will cease to be everywhere regular. However, the gen- 
eralized Painleve-GuUstrand coordinates which are the 
exterior for the non-static cases are everywhere regular, 
and therefore the entire collapse process is everywhere 
well defined. 

By considering only the interior, perfect fluid region, 
and reversing the direction of the time coordinate the 
solution is a cosmological model describing the expansion 
of the universe. 

The article is set-up as follows. In section |n] we de- 
scribe the method for deriving the reduced Einstein field 
equations. Section UTTl analvses the initial and boundary 
conditions required to solve this set of equations. The 
coordinate transformation to put the metric into a gen- 
eralized version of the LTB coordinates is given in IIVI 
and the apparent horizon is briefly discussed in We 
look at the independence of the pressure on the tidal 
forces for the entire spacetime in I VII and finally reduce 
the system to known solutions in section rVIII 

Geometrized units are employed throughout whereby 
c = G = 1. Greek indices run from . . . 3, Latin from 
1 ... 3 and we follow index conventions of Q • 



II. REDUCED EFES 

The metric representing a spherically symmetric space- 
time can be expressed in 3 -I- 1 form, without loss of gen- 
erality, as 
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ds^ = -a^dt^ + Y^iPdt + drf + r'^dn'^, (1) 

where a{t, r) > is the lapse function, P{t, r) is the ra- 
dial component of the shift vector, E{t,r) > —1 and 
dfl"^ := dO"^ + si-a^ ed(j)'^ . 

The energy-momentum tensor for a perfect fluid 
is given by^ 



(2) 



which is related to the Einstein tensor via Gfj,u = SnTu^. 
Here, p is the energy density, P is the pressure and 
is the vector field tangent to the fiuid, which is timelike 
and normalized 

n"n„ = -1. (3) 
Further demanding this satisfies Frobcnius' theorem 

n[^^,,n„] = 0, (4) 



where square brackets denote antisymmetrization, im- 
plies the normal vector is hypersurface forming. The nor- 
mal being timelike implies the hypersurfaces are space- 
like. A particular solution to equation Q is 



(5) 



where t is the temporal coordinate which implies the nor- 
mal one-form can be written in component form as 



-a, 0,0,0] . 



(6) 



The 3 + 1 formalism allows for great simplification 
in spherically symmetric spacetimes. In particular, 
consider any trace-free, symmetric, spatial two-form, 
Wij say. Component-wise, spherical symmetry implies 
We^ = W^'^ . Furthermore, the tensor being spatial im- 
plies Wt = 0. Now, as the tensor is trace-free, the radial 
component must be related to the angular components 
according to Wg^ ~ —2WJ. Therefore, one can write 

= w(t,r)P„-, (7) 

where w{t,r) is the distinct eigenvalue of Wij, and 



Pr 




(8) 



Many of the 3-1-1 variables can be expressed in terms 
of their eigenvalues, analogously with equations and 



• The trace-free extrinsic curvature 



Aij := Kij - - _Ly K := a{t,r)Pij. 



• The trace-free three-Riemann tensor 

1 
3 



(9) 



=^Q,^- := ^R,j - ^ := q{t, r)P,,. (10) 



• The trace-free Hessian of the lapse function 
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1 



D,Dja - — D'Dka := e(t, r)P,j 
6a 



(11) 



Here, Kij and K are respectively the extrinsic curvature 
and its trace, J-ij is the three-metric on spatial hypersur- 
faces, '^Rij and ^R are the three-Ricci tensor and scalar 
respectively, and Di is the unique metric connection as- 
sociated with _Ly . 

The conservation of energy-momentum and ADM 
equations can now be expressed as scalar equations rather 
than complicated tensorial relations. The conservation 



^ For more detailed discussion of r elativistic hydrodynamic equa- 
tions, refer the review article, llTl . In particular, one can further 
split p into a rest-mass density component and specific energy 
density component. 



^ The terms defined here are expressed in terms of their metric 
coefficients in the appendix. 
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equations come from the twice contracted Bianchi iden- 
tities which can be expressed in terms of the energy mo- 
mentum tensor 



v„r"„ - 0. 



(12) 



By decomposing this onto and orthogonal to the spacehke 
hypersurfaces, one can derive an equation which relates 
the extrinsic curvature to the comoving derivative of the 
energy density 



Cnp ={p + P) K, 



(13) 



where £„ denotes Lie differentiation with respect to the 
normal vector field. The other equation coming from the 
Bianchi identities is Euler's equation, which reduces to 
a single component relating the pressure gradient to the 
lapse function 



dP 
dr 



jp + P) da 
a dr 



(14) 



An interesting aspect of equations H13|) and p4|l is that 
the equation of state P = —p, implies both p and P are 
constants. This case is equivalent to a vacuum solution 
with cosmological constant, and while this is interesting 
itself we shall exclude it from the remainder of the anal- 
ysis. 

The ADM system consists of ten equations separated 
into four constraints which are satisfied on all spacelike 
hypersurfaces, and six evolution equations. In spherical 
symmetry, these equations reduce to just two constraint 
and two evolution equations. The two constraint equa- 
tions are the Hamiltonian constraint 



In general, Einstein's field equations are second order 
in the metric coefficients. However, it is noted in equation 
(119(1 that all second order derivative terms cancelled one 
another to leave a first order differential equation. While 
this was the case for the momentum constraint, it does 
not occur for the remainder of the equations. However, 
by combining linear combinations of the remaining equa- 
tions, we can search for combinations of terms whereby 
higher order derivative terms vanish. For instance, the 
combination '^R + 12q gives a zeroth order expression in 
E. Therefore, adding equation l|17|) to six times equation 
(|18|l gives an algebraic expression for E in terms of first 
order derivatives of other metric coefficients. 



(9 

= 2r (H- £;) — Ina - SnPr^ 
or 



2r£, 
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(21) 



Substituting this back through the Hamiltonian con- 
straint ((15|l implies the metric coefficients can be related 
to the energy density. This equation contains a single 
radial derivative of the metric coefficients and pressure 
term 



. 2 d 
47rpr = — 
or 



Integration of 
"mass" function 



(1 



d 

E) —hia 
Or 



a 



(22) 



implies the left hand side becomes a 
^ measured at some radius 



(23) 



M{t,r) Att I pa^da. 



+ -K^ - Ga^ = IGnp, 



and the momentum constraint 



d_, 3. _ ~r^ dK 
dr 3 dr 



(15) 



(16) 



After some algebra, the evolution equations can be shown 
to be 

2£nK - -^R - - 9a^ + -D'^Dka =247rP, (17) 
2 a 



£„a — aK + e — q =0. 



(18) 



Now, we would like to utilize the ADM equations to 
express the metric in terms of physical variables. By 
writing out the momentum constraint (|16|l in terms of 
the metric coefficients (see appendix) , a first order differ- 
ential equation on E and a results 



-1 



1 + E 



2(3 da 
o? dr ' 



where operating on a scalar is 

" adt a dr 



(19) 



(20) 



The integration of the right hand side of H22|l requires the 
consideration of the boundary at r = 0. In particular, a 
physical requirement for the treatment of a stellar sys- 
tem is that the radial derivative of the pressure vanish at 
r = 0. Along with Euler's equation and the regularity 
of /3 and i? at r = 0, this implies the first term on the 
right hand side vanishes at r = 0. Regularity of /?, a and 
P at r = further implies that the remaining two terms 
in the right hand side also vanish at r = 0. Therefore, 
the integration process, and the substitution of Euler's 
equation yields 



M 



47rPr = £, 



/? \ + E dP 
'a [p + P) dr^ 



(24) 



where p 7^ —P. 

Substituting equation through equation re- 
sults in an algebraic expression for the lapse, shift and 



* We note that this is not a physical mass quantity, however the 
quantity arising in a natural way suggests to the authors it is of 
some importance. Furthermore, in the dust limit this becomes 
the familiar mass of the LTB dust solutions, and in the vacuum 
limit becomes the Schwarzschild mass. 
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E 



E- 



2M 



(25) 



While the metric being exphcitly dependant on the lapse 
function is necessary due to Euler's equation, equation 
(|25|l implies we no longer require explicit dependence on 
the shift function. Instead, everywhere we substitute the 
shift function /3, for 



(3 = a 



2M 



E, 



(26) 



where we have taken only the positive root. Consistently 
taking the negative root gives an expanding model rather 
than a contracting model. This can be seen as changing 
the sign of the shift vector is equivalent to reversing the 
time coordinate (i.e. t — s- —t). 

Putting equation (|26|l back through equation l|19|) we 
find an evolution equation for the energy function relat- 
ing to the pressure gradient 



-}-C„E^2{^,E^^?f. (27, 

1 + E \ r p + F or 

Substituting (|26|l into equation (|24ll we find 



„ 2M 1 + EdP M 

Cn\ + E=—-— + — +^nPr. (28) 

V r p + F or 

Finally, by expanding out the Lie derivative in the above 
equation and utilizing the evolution equation on the en- 
ergy function, we can reduce this to a much simpler equa- 
tion 



C^M = ATTr^P\l—+E. 

r 



(29) 



Summarily, the system can be expressed as the line 
element along with two evolution equations 



E dt + dr 



l + E 



r^dn^, 

(30) 



dE 
'dt 



I2M ( dE l + EdP\ , , 

aw +£;(_ + 2 -_)=0, (32) 



dr 



p + P dr 



where p ^ — P, and the lapse function satisfies Euler's 
equation 



1 da 

a dr 



-1 dP 



P dr 



(33) 



As it stands the system is underdetermined and must be 
closed by providing an equation of state 



(34) 



Once the equation of state is prescribed, equation ll^ 
can be integrated to find the lapse in terms of the energy 
density, and then the coupled equations H31|l and 132(1 . 
can be solved simultaneously. 

We make a quick note here regarding the remaining 
coordinate freedom associated with the line element given 
in ((3()|l . It is straightforward to see this is invariant under 
the usual rotation of the two-sphere, and a re-scaling of 
the time coordinate. This is pertinent for the following 
sections, as one can see when the radial derivative of the 
pressure vanishes, equation H33|l implies that the lapse 
function is simply a function of the temporal coordinate. 
However, using the freedom to re-scale this coordinate, 
this function can be scaled to unity. 



III. INITIAL AND BOUNDARY CONDITIONS 



Initial Conditions 



Equations I(31I34|I are intrinsically difficult to solve an- 
alytically, in part due to the freedom associated with 
specifying an equation of state. Whether solving analyt- 
ically or numerically, these equations require the input of 
initial and boundary conditions. As we are establishing a 
physical problem, it is desirable to simply prescribe two 
quantities on the initial hypersurface 



p(0,r) 



and -— 



d_p 

dt 



(35) 



Through the definition for the mass (|23() , this implies the 
initial mass and its time rate of change are known. Fur- 
thermore, an equation of state implies the initial pressure 
and its time rate of change are known, and through Eu- 
ler's equation (|33|l . the initial lapse function is known. 
Equation (|31|) can be rearranged such that 



E = 



1 



dM 



a [47rr2 (p -I- P)] dt J 



V 2M 



(36) 



implying the energy function on the initial hypersurface 
is also determined. Thus, prescribing the two quantities 
associated with the density on the initial hypersurface 
provides sufficient initial conditions. 



B. Boundary Conditions 

Before discussing the free boundary condition we first 
look at the simpler interior and exterior conditions. 

At r = 0, all functions are required to be regular for a 
finite time in the evolution. This will cease to be correct 
once the matter has collapsed to a point, at which point 
only the Schwarzschild solution remains. Extra to the 
regularity condition is that the radial derivative of the 
pressure vanishes at r = 0. 

There are many ways to utilize boundary conditions 
at the interface between the matter filled interior and 
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the vacuum exterior. If we were to simply let the en- 
ergy density, as well as the pressure vanish at a single 
free-boundary, then in general, the interior and exterior 
do not continuously match. This is seen by looking at 
equation l|32() and Eulcr's equation H33() . which imply the 
energy function E, and the lapse function a, are respec- 
tively not defined unless the radial derivative of the pres- 
sure vanishes faster than p + P. We note that simple 
equations of state, for example a linear equation of state, 
do not satisfy this condition. We can therefore always 
prescribe this condition as a boundary condition for the 
equation of state. Alternatively, we can define the equa- 
tion of state such that the pressure drops to zero faster 
than the energy density. In this way, there is a finite 
region consisting of dust, such that the pressure has be- 
come negligible. Physically, this represents a sparse at- 
mosphere around the collapsing object and is consistent 
with the concepts introduced in Q to explain phenom- 
ena such as white dwarfs. Of course, one does not have 
to impose this condition, in which case the vacuum will 
not match continuously to the interior. 

The mathematics that has come from providing the 
atmosphere has prompted reasonable physics. The ob- 
server for the spacetime is travelling with some velocity 
relative to the coordinates given by the shift function. 
We can imagine the transition as the observer travels 
from the vacuum region into the matter region. If the 
observer is suddenly confronted with a region of perfect 
fluid with non-zero pressure, then the velocity of this ob- 
server will contain a jump discontinuity. However, if they 
first travel through an atmosphere, then both transitions 
may be made continuously. 

We therefore have two free boundaries for the space- 
time. The first, denoted rQ-^{t), is the interface defined 
by the point at which the pressure vanishes and energy 
density remains. The second free boundary, denoted 
rg^ (t) > rg-^ (t) , is the interface defined by the vanish- 
ing of the energy density. In this way, our equation of 
state is now given by 



That is, 



lim 

-+r-a(t)- 



E{t,r) 



\im E{t,r) 



(38) 



This boundary condition can then be used to deter- 
mine the unique coordinate system for the Schwarzschild 
spacetime that describes both the interior and exterior 
regions in a single coordinate patch. 



IV. GENERALIZED 
LEMAITRE-TOLMAN-BONDI 

General inhomogeneous dust solutions are most com- 
monly expressed in LTB 1, 2, 3] coordinates. We there- 
fore find the coordinate transformation such that the per- 
fect fluid system is in coordinates that generalize the LTB 
model. By letting {t,r,6,<j)) i-^ {T, R,9,(j)), such that 
t = t{T) = T and r r(r, R) where 



dr 
df 



2M 



E 



(39) 



where now a, M and E are all functions of T and R, the 
metric becomes 



1 dR^+r^dn'. (40) 



l + E \ dR 



Furthermore, this coordinate transformation can be put 
through equations H31|l and (|32|l to give 



dM 
'df 



dr dE _ ^l + EdP ^ 
dRdf ~ p + P'dR.' 

and Euler's equation reads 



2M 
r 

-a 



E, 



2M 



E, 



(41) 
(42) 



r / (p,P) for rG (0,raJ 
P for r e [rai,raj 



(37) 



The inner of the two free boundaries rg^ {t) , can now be 
determined by specifying an initial radius rg-^(0), where 
P goes to zero. Demanding p (hence M), and E be 
continuous, and P = 0, in equations H31|l . (|32() and H33() 
enables the evolution of this boundary to be determined. 

The outer free boundary, rgr^{t), is given such that 
p = P = for all r > rg^ it) . This implies the mass 
function is a constant, and furthermore, one can show 
that this mass in the exterior region must necessarily 
be equivalent to the exterior mass on the initial hyper- 
surface. Thus, the mass in the exterior is given by the 
Schwarzschild mass. 

The boundary condition on both free boundaries for 
the energy function is the simple demand of continuity. 



1 da 

'aim 



-1 dP 



P- 



PdR' 



(43) 



While the system in the (i, r) coordinates described the 
black hole region of the spacetime, and the time reverse 
(i.e. t —t) described the white hole region, the (T, R) 
coordinates now describe both these regions. Further- 
more, the Jacobian of transformation is zero only when 
dr/dR vanishes. 

It is apparent how the system of equations H39I43|I re- 
duces to the LTB system as the pressure vanishes. This is 
therefore a direct generalization of the LTB model to in- 
clude inhomogeneous pressures. Furthermore, by looking 
at this reduction of the system, one can appreciate the 
function E, is simply the energy function of the stan- 
dard LTB model. In this way, the cases of E being 
positive, negative or equivalently zero correspond to un- 
bound, bound and marginally bound models respectively. 
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It is interesting to note that the mass function is still 
in terms of the original definition given in equation (|23|l . 
which now reads 



rr{T,B.) 

M{T, R):=Att i pa^da. 

Ja=0 



(44) 



This can be differentiated with respect to the LTB radial 
coordinate _R, and rearranged to give 



p{T,R) 



M' 



(45) 



where a prime denotes differentiation with respect to R. 
This is the standard eq uation for the pressure-free LTB 
solution (see for e.g. and indicates shell-crossing 

singularities occur for r'{T,R) — and shell- focussing 
singularities for r{T,R) — 0. The occurrence of these 
singularities is now implicitly dependant on the pressure 
due to equation 

An interpretation of these singularities is awkward due 
to the choice of coordinates, implying they are exhib- 
ited simply as an infinite density. However, in the {t, r) 
coordinates, the singularities are exhibited as a multi- 
valued mass function. This is akin to fluid shock waves 
(c.f. pressure-free case 

MM) 

Further analysis of these 
shocks will appear in a future article. 



V. APPARENT HORIZON 

To establish the apparent horizon for a spacetime, one 
must analyse both the metric as well as the equations of 
motion. In particular, utilizing the unique afBne parame- 
ter associated with the metric. A, and the Euler-Lagrange 
equations, there are five equations governing the dynam- 
ics of all null geodesies. 



^ _ -a^l-2M/r) .^ ^ 2ay/2M/^ 



E 



l + E 



„2/)2 



= 



^ l + E ' ' " 
d dC _ dC_ 
dX dif^ dxt^ ' 



r'^ sm 



l + E 
2 0<i)'^ = 0, 



-tf 



(46) 
(47) 



where a dot represents differentiation with respect to A, 
and C — C{x^,x^) is the Lagrangian. Spherical sym- 
metry implies only the radial geodesies are required, and 
we therefore set 9 = cf) — Q. By utilizing the remain- 
ing three equations, one can solve for t and f with one 
arbitrary constant of integration. We next define a null 
vector according to 



Ilx'' 



(48) 



which is therefore everywhere tangential to the congru- 
ence of radial null geodesies. The arbitrary constant acts 



to scale this vector, and it is therefore set to unity with- 
out loss of generality. We find 



A:^ = ^/lTE 



-,VT+E 



2M 



E,0,0 



(49) 



where the radial coordinate is now parametrized by the 
temporal coordinate. It is straightforward to verify this 
vector is null, and the expansion of this vector field is 
simply it's divergence. 



(50) 



This is a measure of the convergence and divergence of 
the congruence of radial, null geodesies. The limiting 
case of the converging and diverging geodesies are those 
given by a vanishing expansion factor, = 0, and this 
defines the apparent horizon l22j. Evaluating this term 
gives 
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VTTE d 
OLT^ dr 



ar 



2 VT+E- 



2M 



E 



(51) 



By setting = 0, integrating from r = to some finite 
radius and using the regularity of the functions at r = 
implies the apparent horizon is given by the parametric 
equation 



r(i) = 2Af (i,r(t)). 



(52) 



This simple form shows that at the interface between 
the matter and vacuum regions, where the mass func- 
tion simply becomes the Schwarzschild mass, the appar- 
ent horizon reduces to the familiar event horizon in the 
Schwarzschild spacetime. 



VI. TIDAL FORCES 

Ellis [IJl showed that the tidal forces of a spherically 
symmetric spacetime are given by the electric confor- 
mal curvature Eij, which arises from the decomposi- 
tion of the Weyl tensor. The electric conformal cur- 
vature is a trace-free, spatial two-tensor, and therefore 
has a unique eigenvalue, denoted A. Spherical symme- 
try further implies this is the only non-zero contribution 
from the Weyl tensor. A constraint equation from the 
gravito-electromagnetic system of equations relates the 
divergence of the electric curvature to the gradient of the 
energy density 



(53) 



Evaluating this in terms of its eigenvalue implies the 
equation can be integrated, and we find the relation be- 
tween A and p is given by 



A = 



— r d 
3 dr 



(54) 
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It is interesting to note the tidal forces on the space- 
time are exphcitly independent of the pressure (this is 
the same result achieved in Furthermore, the ex- 

terior Schwarzschild spacetime is described by M = Ms, 
and the usual tidal force for the Schwarzschild spacetime 
results. 

While in general these are Petrov type D solutions, the 
system can be reduced to a Petrov type solution. These 
are found when the above eigenvalue vanishes, which in 
turn gives an FRW interior solution. The FRW solu- 
tions have homogeneous matter distributions, implying 
the pressure and density are constant along each t equals 
constant slice. This homogeneity, along with the spher- 
ical symmetry, in turn implies there are no tidal forces 
felt by any particles in the spacetime. 



VII. ANALYTIC SOLUTIONS 

There are many methods available to us in searching 
for analytic solutions of equations (jH 1133(1 . For instance, 
one can employ an equation of state, solve H33|l and at- 
tempt to solve the remaining two coupled equations. This 
process, while its motivations are purely physical, is over- 
run with mathematical difficulties. For instance, choos- 
ing an arbitrary equation of state will, in general, imply 
the time variation equation on the mass is highly non- 
linear. 

We therefore make assumptions which are more di- 
rected at simplifying the mathematics, in order to grasp 
a better handling of the physics. Once the mathemat- 
ics is simplified, we can analyse the resulting equation of 
state to ascertain its physical relevance. 

We have already determined that utilizing an equation 
of state such that P = implies the system reduces to the 
LTB dust interior with generalized Painleve-GuUstrand 
exterior region [l5l |. Another major class of important 
solutions are the static solutions. 



A. Static Solutions 

We retrieve static solutions of the field equations by 
removing all time dependance. From equation 1(31(1 we 
either derive P = —p, which is a static (anti) de-Sitter 
model, or 



E 



2M 



(55) 



This is a statement implying that the shift vector (3, is 
identically zero. The metric is therefore given by 



1 - 2M/r 



+ r^dn^, (56) 



and furthermore, equation ((27() implies E = E{r), which 
implies M = M{r), and in turn p — p{r). Finally, 



putting this back through equation 1(28(1 and utilizing Eu- 
ler's equation 



-1 dP 



AI + iirPr^ 



p + P dr r2 (1 - 2M/r) 



a dr 



which is exactly the TOV equation of hydrostatic equi- 
librium. 

The static solutions are a special case of the equa- 
tions we have derived, and all other cases (that is, with 
E ^ ~ 2M/r), result in non-static systems. While this is 
not a new result, it is interesting to note that letting the 
pressure and density vanish at some finite radius implies 
an exterior region is simply given by the Schwarzschild 
spacetime in Schwarzschild coordinates. To see this we 
let P ~* and p — > (by utilizing the atmosphere dis- 
cussed in section Ullll . and find the solution for the lapse 
function is 



1 



2M, 



(58) 



Putting this into l(56() gives the Schwarzschild metric in 
Schwarzschild coordinates. 

Therefore, the natural coordinates for the exterior 
of the static model are the Schwarzschild coordinates. 
While the generalized Painleve-GuUstrand coordinates 
are everywhere regular, the Schwarzschild coordinates 
have a coordinate singularity at a finite radius, r = 2Ms. 
However, the interior for the static case must necessar- 
ily have rg > 2Ms, and therefore the coordinates for the 
entire spacetime are still everywhere regular. 



B. Non-Static Solutions 

1. Linear Mass Equation 

The form of the mass equation l(31() is interesting. Es- 
sentially, once an equation of state is prescribed, both the 
pressure and the lapse function can be written in terms of 
radial derivatives of the mass function. We therefore ask 
the question: under what conditions will equation 1(32(1 
be a quasi- linear differential equation? That is, with a 
linear radial derivative. This question can be expressed 
mathematically as 



a 



dM 
dr 



dM 
dr ' 



(59) 



where fc is a constant of proportionality. Rearranging, 
and using the definition of the mass function implies 



kp 



(60) 



P + P 

Substituting this through Euler's equation 133() implies 
P dp 



p{p + P) dr 



0. 



(61) 
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Therefore, equation l|31|l is linear either in the absence 
of pressure (see ) or if the energy density is indepen- 
dent of the radial coordinate. Considering a reasonable 
equation of state with this latter condition implies P is 
also independent of the radius. These two conditions are 
exactly an FRW interior with vacuum exterior. 

The right hand side of Euler's equation (|33|l now van- 
ishes, implying the lapse is simply a function of the time 
coordinate. Residual coordinate freedom can be utilized 
to re-scale the time coordinate such that a = 1 without 
loss of generality. This implies only a linear, barotropic 
equation of state is valid 



P=(fc-l)p. 



(62) 



The differential equation on the mass function H31|) can 
now be expressed in terms of the density and energy func- 
tions 



dp 3 /Stt „ ^ , ^, 



(63) 



Functional dependance for this equation is consistent if 
and only if 



(64) 



These solutions are actually characterized by another as- 
pect of the solution. Namely, they form a subset of the 
shear-free solutions. Shear-free solutions are those such 



that Ai 



0, which in the formalism given above is 



equivalent to a{t,r) = 0. In general, these solutions can 
be shown to satisfy y^2M/r + £^ = rf{t), for some func- 
tion of integration f{t). The set of shear-free solutions 
we discuss differ from the shear-free solutions in .6j as 
Aij is a coordinate dependant quantity. 

Returning to the smaller class of FRW solutions, x ^nd 
p now satisfy a coupled system of first order ordinary 
differential equations 



|=3^fp + x(P + P), 



(65) 
(66) 



This system can be solved for the interior region using 
the equation of state (|62(l . to give 



(67) 



where A is an integration constant. This system can fur- 
ther be solved analytically when the constant fc, is spec- 
ified. Furthermore, at some finite radius on the initial 
hypersurface, the density and pressure will be zero. In 
this region, equation (|66|l can be calculated to find the 
unique exterior region of the spacetime that matches onto 
the FRW interior. 



2. Spatially Quasi-Flat Solutions 

We can also search for solutions that contain some 
sense of spatial flatness. Letting the three-Ricci tensor 
vanish implies spacelike hypersurfaces are flat. This con- 
dition implies the energy function necessarily vanishes, 
which in turns gives P = P{t). This is therefore equiva- 
lent to the FRW solutions discussed above. 

Alternatively, we can let the three-Ricci scalar '^R, van- 
ish. This implies E — f{t)/r, where f{t) is a function 
of integration. Without loss of generality, we can utilize 
residual coordinate freedom to re-scale the time coordi- 
nate such that / = 1. The energy function is therefore 
given by 



r 



(68) 



Therefore, the energy plays a similar role to the mass 
function as they are both represented via terms which 
are inversely proportional to the radial coordinate. By 
rewriting equation (|32|l . and using Euler's equation, we 
find 



1 



1 



2r (r + 1) p + P dr 



dP -1 da 
a dr 



= -in- (69) 



Without an equation of state, this equation can be inte- 
grated for the lapse function, yielding 



a = g(i)\/l + J, 



(70) 



where g{t) is a bmction of integration. Furthermore, by 
just specifying P = P{p), we find a formal integral for 
the energy density 



dP dp 



P {p) + p dp dr 



dr = In 



1 



(71) 



Upon specifying an equation of state, equations (|70() and 
()71|l can be substituted back through the evolution equa- 
tion for the mass H31|) to establish constraints for the 
arbitrary function g{t). 



3. Self-Similar Solutions 

While self-similar fluid solutions have been readily 
treated in the literature (for a recent review see Q), it is 
a worthwhile exercise to establish how they arrive in the 
context of the equations derived herein. 

Spherically symmetric self-similar solutions can be put 
into a form where all dimensionless quantities are func- 
tions of the self-similar variable ^ := t/r For instance, 
we define H {£,) := 2M/r and p{^) := AirPr"^. Further- 
more, a, E and /3 are all dimensionless, and are therefore 
functions of the self-similar variable ^. This implies the 
density becomes 



d^ 



(72) 
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Equations H31|l . (|32|l and (|33() respectively become 
dH P{H + 2p) 



1 dE 



-4/3 



1 + E d^ ^i{l + |30 \d^ 
1 da —2 / do 



where 



and 



d^ 



P (0 = aVH + E. 



(73) 
(74) 
(75) 

(76) 
(77) 



The system has therefore reduced to a comphcated set 
of coupled ordinary differential equations. Although the 
system is now written as ordinary differential equations, 
it is apparent that they are still highly coupled, and 
therefore difficult to solve analytically. However, physical 
aspects of the solutions can be explored through methods 
of dynamical systems. 



VIII. CONCLUSION 

We have derived reduced field equations describing a 
spherically symmetric spacetime with a perfect fluid re- 
gion in the interior and an exterior vacuum region. The 
interior is shown to be a direct generalization of the LTB 
dust solutions to include inhomogeneous pressures. A 
general solution to these reduced field equations is not 
found, and it is not obvious that one such solution exists. 
Despite this, we have shown a handful of specific solutions 
that were already known. It is part of the robustness 
of the formalism that everything from the static, TOV 
equations, as well as the FRW models, the self-similar 
solutions and many more known solutions fall naturally 
from the equations derived herein. 

Before contemplating exact solutions of the reduced 
equations, an equation of state must be prescribed. 
Within the present work we have made no attempt to 
prescribe such an equation due to the unreliability of the 
present understanding of microscopic physics under ex- 
treme pressures and densities. Rather than prescribe an 
equation of state, it is a common technique to use ge- 
ometrical or mathematical simplifications such that the 
field equations are put into a form that is soluble (see 
for a review). From this point, mathematically plau- 
sible equations of state are derived. While this method 
yields tractable solutions, it is still an open question as 
to whether the derived equations of state are physically 
reasonable for high density and pressure regimes. 

The above method is often used to analyse the TOV 
equation describing the static case. This is because only 
a handful of analytic solutions are known with simple 



polytropic equations of state. While we claim no new 
ground on this solution, we point out that it arises natu- 
rally from our formalism. The TOV equation is the only 
case that has a diagonal metric, implying it is also the 
only case where the Schwarzschild coordinates provide a 
natural description for the exterior region. 

An interesting aspect of this is that the solution is al- 
ways regular through the horizon. This is obvious in the 
collapsing cases, whereby both the interior and exterior 
regions are described by coordinate systems that are ev- 
erywhere regular (r ^ 0). However, we showed the static 
case reduces to the diagonal metric, with Schwarzschild 
coordinates in the exterior. This system is still every- 
where regular, as a static object must necessarily have 
fluid covering r — 2Ms. 

There are many avenues for future work within the 
realms of this formalism. Obviously, solving the system 
of derived differential equations both analytically and nu- 
merically with prescribed equations of state is a desirable 
research goal. Another aim is to study the formation and 
evolution of both shell-crossing and shell-focussing singu- 
larities. While this has been extensively researched for 
the dust models, very little is understood about the cases 
with the inclusion of pressures. 

The model we have described is also able to be gen- 
eralized through many avenues. Firstly, the inclusion of 
heat flux and anisotropic pressure terms will allow for 
more realistic stellar models. With only a perfect fluid 
interior and Schwarzschild exterior, the model is akin to 
a pressure cooker. As the system collapses, the temper- 
ature will naturally rise, however, there is no allowance 
within the model for the inclusion of radiation in the 
form of either photons or gravitational waves. The inclu- 
sion of aniostropic terms for the interior, and generalizing 
the exterior to the Vaidya spacetime will allow for inco- 
herent null radiation (see for e.g. [S^, [l^l, and 
references therein). Relaxing the spherical symmetry to 
quasi-spherical symmetry will allow for gravitational ra- 
diation where the exterior will be a Robinson- Trautman 
spacetime. One final achievable generalization of the 
model presented herein is to derive fluid equations for 
the axisymmetric case. This will include a rotating fluid 
for the interior with a Kerr exterior. 
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APPENDIX 

Terms in the ADM equations are expressed as func- 
tions of the metric coefficients for clarity 

'R=^iirE), (A3) 



^^-rVTTE^d_/y/TTE_da\ 
3a dr \ r dr J ' 

'-D>^D,a=:^^(r^VTTE'^]. (A6) 
a ar^ or \ or I 
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